A toy model of faith-based systems evolution 
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A simple agent-based model (ABM) of the evolution of faith-based systems (FBS) in human social networks 
is presented. In the model, each agent subscribes to a single FBS, and may be converted to share a different 
agent's FBS during social interactions. FBSs and agents each possess heritable quantitative traits that affect the 
probability of transmission of FBSs. The influence of social network conditions on the intermediate and final 
macroscopic states is examined. 
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I. INTRODUCTION 

Human societies tend to hold sets of shared and interrelated 
beliefs that are not directly or easily accessible to proof or dis- 
proof. Such sets of beliefs, which will be called Faith-Based 
Systems (FBSs), may include, for example, cultural prefer- 
ences, personal worldview, and belief or disbelief in one or 
more supreme beings. FBSs are often highly detailed and tend 
to be shared within local communities, but are highly diverse 
among global communities. 

Dawkins |1] suggested that the evolution and spread of 
FBSs and other cultural information (memes) may be subject 
to natural selection, noting that both parasites and memes are 
transmissible between hosts and show heritable variation, cre- 
ating the potential for competition for survival and reproduc- 
tion. While the meme hypothesis has received considerable 
theoretical attention, there have been few attempts to explic- 
itly model it. 

Previous models of social transmission j^|^|^|5|] have ex- 
amined mainly two-state systems; in such models, agent opin- 
ions are considered to be transmissible Boolean states. Such 
Boolean state models do not allow for the evolutionary poten- 
tial of memetic transmission, which may influence final state 
distribution, or for the potential for a large number of FBSs 
to exist simultaneously in a single social system. In contrast, 
FBSs can be modelled as continuous states on a binary vector. 
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allowing for memetic evolution among populations of varying 
FBSs. 

In this paper an agent-based model of the evolution and 
spread of FBSs in a two-dimensional lattice is presented. 
Agent-based models are useful for examining how local prop- 
erties of individual units interact to produce macroscopic 
properties of systems, and have been applied in the analysis 
of social opinion formation cultural interactions 

|6], epidemiology and population dynamics ff\. The current 
model explores the development of FBSs with evolutionar- 
ily relevant properties (faith, proselytism and mutability), in 
simple, individually varying agents. Mortality and fecundity 
probability curves for agents are applied based on human pop- 
ulation data from the 1990s |8]. Agents have heritable per- 
sonality characteristics that affect their susceptibility to con- 
version and their likelihood of creating a new FBS; this pro- 
vides a non-uniform spatial environment for the evolution of 
FBSs. Characteristics of agents modelled are (cf. Tab.UJi re- 
sistance (R) (the tendency to resist conversion to a different 
FBS), activity (A) (the tendency to participate in transmission 
of FBSs), charisma (C) (ability to convert other agents to their 
current FBS), and imagination (/) (the tendency to invent a 
new FBS based on the existing one). 

The influence of two factors on the distribution and charac- 
teristics of FBSs is examined. The frequency of social inter- 
actions of agents, where social interactions offer the potential 
for conversion between beliefs, seems likely to influence the 
spread and survival of FBSs. When social interactions are less 
frequent, pressure for conversion between FBSs is likely to be 
smaller, but social support for existing FBSs is also decreased. 
If social support plays an important role in preventing conver- 
sion, extinction of unsuccessful FBSs may occur more quickly 
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when social interactions are less frequent. This could lead to 
decreased numbers of FBSs, each on average having a larger 
number of adherents. Additionally, the effect of varying total 
social network size on FBS evolution is explored. Larger total 
social network size may increase the influence of the social 
interaction condition, and may allow a larger number of coex- 
isting FBSs and older FBSs (due to lower extinction rates in 
larger grids). 

II. MODEL 

The model initiafly consists of iV = 10 x 10 or iV = 20 x 20 
agents which are located on a regular rectangular two dimen- 
sional lattice with free boundary conditions. Each agent, oc- 
cupying a single site on the grid, interacts with up to eight 
social neighbors which are randomly chosen from nearest ge- 
ometric neighbors; agents located at the border of the grid 
have fewer neighbors, (cf. agent k on Fig.QJ. For comments 
on other possible grids see 
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FIG. 1: Sample agent m with its 3 x 3 environment. 1, 4 or 8 in- 
teracting agents (social neighbors) are randomly chosen from 3x3 
geometrical environment. 

Within the model, social space is created by a regular rect- 
angular lattice and social distance is equal to one for the near- 
est interacting neighbors and infinity for others, i.e. there are 
no long range interactions (in strict meaning). This may be 
somewhat analogous to social transmission in pre-literate hu- 
man societies, but is less applicable to modern societies. 



AGENT 



Believe 


B 


FBS accepted by a given agent 


Activity 


A 


Probability of participating in social interactions 


Charisma 


C 


Ability to convert other agents 


Imagination 


I 


Ability to create a new FBS 


Resistance 


R 


Resistance to conversion by other agents 



TABLE I: Characteristics of agents. 

Every agent is characterized by a set of internal parameters 
which are constant over the life time of an agent (cf. Tab. 
Features of agents provide a non uniform environment for the 
evolution and spread of FBSs. 



Every FBS is characterized by a set of external parameters 
which are constant over time (cf. Tab.lIH. Since a large num- 
ber of agents may share the same FBS, a single parameter 
set can be widespread and sometimes dominant throughout a 
grid. 



FAITH BASED SYSTEM (FBS) 



Faith 


F 


Resistance to conversion pressure from 






other FBSs 


Proselytism 


P 


Strength of conversion pressure against 






other FBSs 


Mutability 


M 


Probability of mutation of a given FBS 



TABLE II: Characteristics of FBSs. 

All characteristic parameters of FBSs and individuals are 
randomly generated from a restricted Gaussian distribution 
with mean /i = 10 and standard deviation ct = 3 and con- 
verted into integers. The above restrictions are related only 
to the initial state distribution; over time values of parameters 
can change due to mutation processes. 

The parameters in the model (Tab.lllandlllli reflect a limited 
selection of features of real agents and FBSs. These parame- 
ters were chosen to express some of the most influential prop- 
erties of real agents and FBSs, while limiting the parameter 
space explored. 

Agent TO receives support S from social neighbors (sur- 
rounding interacting agents) (cf. Fig.Q which accept the same 
FBS B{m) 

S{B{m)) = SBi,n)Ms)G{s)P{B{m)). 

s^{i ,j ,k ,1 ,n,p ,q ,r} 

(1) 

where 6ij is a Kronecker delta function, P{B{rn)) is pros- 
elytism of the FBS accepted by the agent m and all other 
symbols are described in Tab. |I] and HH The agent to is also 
subjected to the environmental pressure V (cf. Fig.Q) which 
is caused by surrounding interacting individuals having differ- 
ent FBS, e.g. pressure caused by agent k accepting FBS B{k) 
{B{k) ^ B{m)) 

V{B{k))= 5Bik),Bis)C{s)P{B(k)). (2) 

The sum in eq. ([0 and (|2} is performed over s G 
{i,j,k,l,n,p,q,r} i.e. s takes 1, 4 or all possible values 
from the set {i, j, k,l,n, p, q, r}. The Sij in eq. Q and (|2j 
with appropriate indices allow to choose agents sharing the 
same FBS, which in particular may be different than accepted 
by the agent to. Support S and pressure V (eq.[2and|3 are 
weighted sums of social neighbors' charisma as a weight the 
proselytism of FBSs are taken. 

Every agent is characterized by total self resistance TZ 

n{m) ^ R{m)F{B{m)). (3) 

Total self resistance is calculated as a product of an agent's 
resistance R{nn) and faith of an agent's FBS F{B{m)). 
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Social interactions between agents cause the agent m ac- 
cepting FBS B{m) to be converted into the FBS B{k) of an 
agent k if 



V{B{k)) > max[S{B{m)), TZ{m)] 



(4) 



Transmission of FBSs can be interpreted in terms of epidemi- 
ology as infection, in that an FBS depends on a host agent, 
and can be transmitted between agents via specific forms of 
contact. In our model, transmission is limited both by the in- 
dividual resistance of an agent, and by the presence of social 
support for the agent's current FBS. This attempts to reflect 
the interaction of personality and environment in social trans- 
mission; for some individuals, personal characteristics are the 
main factors preventing conversion, while for others, social 
support is more important. 

New FBSs are created randomly with probability dependent 
on the imagination of the agent and the mutability of its cur- 
rent FBS, e.g. agent m can create new FBS with probability 
oc I{m)M{B{m)). The new FBS has the same characteris- 
tics as the old FBS, with a random increase or decrease by 1 
unit to any one of its characteristic values (proselytism, faith 
or mutability). This mechanism introduces an evolutionary 
principle, which allows directed change in FBS characteris- 
tics over time. 

The probability that a given agent (e.g. m) is going to par- 
ticipate in social interactions during a given iteration of the 
model is oc A{m). Additionally, during each iteration, agents 
can move to a randomly selected adjacent cell, provided it is 
empty. 

Space occupied by agents is limited to the total area of the 
grid, but the number of agents can vary over time due to mor- 
tality and fecundity of individuals. The agents can reproduce, 
producing a new agent with similar personality characteristics 
in any adjacent empty cell. The probability of reproduction 
depends on the age of the parent agent and is proportional to 
121 - (< - 25)2 with t e [14, 36]. The agents die at ran- 
dom according to a probability curve based on data for Aus- 
tralian, English, Polish and US populations during the 1990s 
fs]. Conditional probability to die at given age t is propor- 
tional to {t — 12)2. xhese phenomenological relations reflect 
the shape of real life fecundity and mortality curves. More 
details about the phenomenological theory of mortality and 
aging can be found in lIlOll . 

In the model the agents are updated in a random order. 
This kind of updating corresponds to local social interactions 
which, as in reality, occur non-synchronously within small so- 
cial groups. Social interactions are governed by a determinis- 
tic rule applied to randomly chosen agents. These dynamics 
create a modified version of social impact theory flT'l . 



III. RESULTS 

Macroscopic features of the intermediate and final state for 
different lattice sizes (10 x 10 and 20 x 20) and different num- 
bers of interactions per agent per timestep ( 1 , 4, 8 agent inter- 
actions) (cf. Fig.Q were analyzed. Data were collected every 
1000 iterations until the final state was reached (state after 



10000 iterations). To obtain sufficient statistics, around 1000 
simulations were performed for each model considered. Prob- 
ability density functions (PDFs) of number of FBSs (Fig.|2jl, 
age of FBSs (Fig.0} and numbers of adherents (size of FBSs) 
(Fig-EJ are shown. Additionally, time behavior of these vari- 
ables (Fig. |3and|5} and correlations between FBSs' charac- 
teristics were examined. 

PDFs of a number of FBSs seem to have exponential tails 
(cf. semi-log Fig.|2j, i-e. 



P(#FBS) oc exp(-a#FBS). 



(5) 



There is a high probability of a small number of FBSs, and a 
small probability of a large number of FBSs co-existing in a 
single system. Values of exponents a for the lattice 10 x 10 
are higher than for the grid 20 x 20. Time behavior of expo- 
nents a (with error bars) is presented in Fig.|3l For t > 4000 
there is a difference between values of exponents a for binary, 
1 to 1, interaction (higher exponent) and group, 4 or 8 to 1, in- 
teraction (lower exponent), for both lattice sizes. Values of 
exponents a for 4 or 8 interacting agents are similar. 
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FIG. 2: Probability distributions (semi-logarithmic scale) of number 
of FBSs at t = 10000, grid size 10 x 10 (left panel) and 20 x 20 
(right panel), for a given number of interacting agents. 



At any point in time it is possible to find an FBS as old as 
the system. There is also a significant probability of finding 
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FIG. 3: Time behavior of characteristic exponents a (cf. eq. |5} for 
number of FBSs distributions, grid size 10 x 10 (left panel) and 20 x 
20 (right panel), for a given number of interacting agents. 



FIG. 4: Probability distributions (logarithmic scale) of age of FBSs 
at t = 10000, grid size 10 x 10 (left panel) and 20 x 20 (right panel), 
for a given number of interacting agents. 



young and very young FBSs. The age distribution expresses a 
power law (cf. log-log plot Fig.0}, i.e. 

P(age) oc age"''. (6) 

Small deviations from a straight line, which are visible on the 
log-log plot are caused by the finite world size and limited 
size of random sample. The order of magnitude is affected 
by the age of the model. Time behavior of exponents f3 (with 
error bars) is presented in Fig. |5] Values of exponents f3 for 
every type of interaction depend on the size of the lattice, but 
do not change dramatically over time (cf. Fig.|5}. For binary 
interactions, values of exponents /3 are smaller than for group 
interactions. Interactions with 4 or 8 agents seem to be equally 
efficient, i.e. values of exponents are similar 

Each agent occupies a single cell on the lattice, i.e. the size 
of FBS is the number of adherents. Distributions of size are 
shown in the Fig.|6] In the log-log plot the linear decay and 
local maxima of PDFs for both lattice sizes are visible. For 
the 20 X 20 grid the linear part of the PDF is steeper than for 
the 10 X 10 grid. For the bigger world size, differences in 
slope for binary and group interactions are observable. The 
differences for the small world size are not so apparent. Lo- 



cal maxima of PDFs which are visible (cf. Fig.|6jl correspond 
to the situation when a large cluster of a single FBS occupies 
almost the whole lattice. The size of a macroscopic cluster is 
approximately 70 — 80% of the whole lattice. Local maxima 
of PDFs are most visible for the 20 x 20 lattice with 4 or 8 
agents interaction and are not visible for binary interactions 
on the same lattice. This phenomenon can be explained by a 
kind of "non-local" character of a social interaction with 4 or 8 
social neighbors. Despite the fact that in the interaction only 
some of nearest geometrical neighbors participate, the inter- 
acting agent is more linked with its environment. The larger 
number of links provides a kind of social connection, which 
introduces effects similar to long range interactions. This is in 
agreement with the theory of social distance, which in general 
can be different from geometrical distance. For lattice 10 x 10 
the local maxima of FBSs' size distribution are visible for all 
interaction conditions. This is probably a consequence of the 
finite world size. 

Effects of the number of interactions per agent are visible, 
manifested as different values of exponents for systems with 
only 1 agent interaction compared with systems with 4 or 8 
agent interactions (Figs. |2] |3] and |3 15}. Binary interactions 
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FIG. 5: Time behavior of characteristic exponents /3 (cf. eq. |6} for 
age of FBSs distributions, grid size 10 x 10 (left panel) and 20 x 20 
(right panel), for a given number of interacting agents. 



FIG. 6: Probability distributions (logarithmic scale) of number of 
adherents of FBSs att = 10000, grid size 10 x 10 (left panel) and 
20 X 20 (right panel), for a given number of interacting agents. 



seem to be more efficient in terms of FBS transmission. This 
efficiency is manifested by a smaller number of FBSs exist- 
ing in the 1 to 1 interaction condition, in comparison with the 
group interaction conditions. Additionally, FBSs in the 1 in- 
teraction condition are older and have more adherents. The 
phenomenon of higher efficiency can be explained by the lack 
of support for agents participating in binary interactions. Bi- 
nary interactions might be seen as corresponding to informed 
discussions in real social interactions, which provide an effi- 
cient way of convincing people by largely avoiding negotia- 
tion of opinions. 

Initially all characteristics were generated as independent 
and identically distributed random variables, according to the 
restricted Gaussian distribution. Over time the mechanisms 
of conversion and mutation included in the model introduce 
coiTelations. Positive correlations (0.58-0.80) between the 
age of an FBS and the number of adherents can be observed. 
Within the model, it seems that even highly successful FBSs 
need time to become widespread. The variations of coiTela- 
tion value are caused by mortality and conversion processes. 
The absolute values of other correlations are smaller In par- 
ticular, small negative correlations (-0.17 — 0.08) between 



mutability and the age of FBSs, and mutability and the size 
of FBS were observed. As mutability increased, both the age 
and size of FBSs decreased; this is attributed to the increased 
tendency for agents holding high-mutability FBSs to invent 
new FBSs. Smaller negative correlations also occurred be- 
tween proselytism and age, and proselytism and size. Because 
of the generally low level of correlation, it is hard to draw 
any detailed conclusion from these findings. What is intrigu- 
ing, however, is an unexpected sign of correlation between age 
and proselytism, and size and proselytism. This effect may be 
caused by the limited size of our random samples or the short 
time of simulation. It is also necessary to emphasize that if 
random variables are independent they are uncoiTelated. The 
inverse theorem is not necessarily true. 



IV. SUMMARY 

A simple agent-based model of the formation and evolu- 
tion of FBSs has been presented. The influence of the lattice 
size and frequency of social encounters on macroscopic out- 
put have been examined in detail. 



6 



In the model, the age distributions of FBSs express power 
laws (Fig. 0}. Larger social networks, as reflected in larger 
lattices, have larger numbers of coexisting FBSs; on aver- 
age, FBSs in such systems are younger and have more adher- 
ents. At low levels of social encounters, FBSs are, in general, 
older and have more adherents, than at higher levels of so- 
cial encounters. A possible explanation of this pattern is that 
the single-encounter condition prevented agents from receiv- 
ing support for their current FBS from adjacent interacting 
agents of the same FBS; this could cause a higher incidence 
of conversion leading to more rapid extinction of unsuccessful 
FBSs. 

On the distributions of the FBSs' size (cf. Fig.|6|l local max- 
ima of PDFs are visible. These maxima manifest the state in 
which almost all agents share the same FBS. These phenom- 
ena are clearly visible for non-binary interactions which in 
social systems correspond to long range interactions. 

The toy model which was presented above allows for fur- 
ther possible extensions. In particular, it would be interesting 
to examine fitness effects of FBSs on agents, different dis- 



tributions of agent and FBS traits, and more realistic social 
network and agent traits. Such extensions may enable a better 
critique of the utility and robustness of models that perceive 
FBSs as objects of natural selection. 
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